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^^ , Abstract 

This paper is concerned with the problem of a semi-infinite crack steadily propagating in an elastic solid 

j___, with microstructures subject to antiplane loading applied on the crack surfaces. The loading is moving 

(~ i ' with the same constant velocity as that of the crack tip. We assume subsonic regime, that is the crack 

CL(' velocity is smaller than the shear wave velocity. The material behaviour is described by the indeterminate 

i . theory of couple stress elasticity developed by Koiter. This constitutive model includes the characteristic 

"^j ' lengths in bending and torsion and thus it is able to account for the underlying microstructure of the 

Cu \ material as well as for the strong size effects arising at small scales and observed when the representative 

S . scale of the deformation field becomes comparable with the length scale of the microstructure, such as 

I I ' the grain size in a polycrystalline or granular aggregate. 

The present analysis confirms and extends earlier results on the static case by including the effects 

\^ ' of crack velocity and rotational inertia. By adopting the criterion of maximum total shear stress, we 

J^ ' discuss the effects of microstructural parameters on the stability of crack propagation. 

C"*~- I Keywords: Couple stress elasticity; Dynamic Fracture; Steady-state propagation; Microstructures; Integral 

O^ ■ transform 

cn . 

O ■ 1 Introduction 

(N; 

^~~^ , In view of the discontinuous nature of many engineering materials, it is well known that the classical theory 

^ ' of elasticity is not adequate to accurately describe the stress and displacement fields at small scales (Morozov, 

J^ , 1984). The discrepancy between the classical theoretical predictions and experimental results is found more 

rS ' pronounced for materials with a coarse-grain structure (Fleck et al., 1994). The mechanical behaviour of 

C^ , most materials with microstructure, like composites, cellular materials, foams, masonry, bone tissues, glassy 

and semicrystalline polymers, is strongly influenced by the microstructural characteristic lengths, especially 

in the presence of large stress (or strain) gradients (Lakes, 1986; Lakes, 1995). These findings stimulated the 

development of generalized theories of continuum mechanics such as micropolar elasticity (Cosserat, 1909), 

indeterminate couple stress elasticity (Koiter, 1964) and more recently strain gradient theories (Aifantis, 

2011; Fleck & Hutchinson, 2001). An extensive review of various formats of gradient elasticity and their 

performance in static and dynamic applications can be found in Askes & Aifantis (2011). 

Because of the characteristic lengths coming into play, the generalized theories predict dispersive wave 
propagation in microstructured media (Nowacki, 1985). Also the generalized theories predict different asymp- 
totic behaviour for some conrponents of the deformatioir and stress fields near singular points, compared to 
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the classical elasticity theory (Morozov, 1984). Correspondingly, the notion of energy release rate and J- 
integral needs to be generalized (Lubarda and Markenscoff, 2000; Piccolroaz et al., 2012). Interestingly, 
the relation between stress intensity factors for classic and CS elasticity discussed in detail in Radi (2008) 
by means of the Wiener-Hopf technique was also found by Mozorov & Nazarov (1980) and by Nazarov &: 
Semenov (1980). The fact that microrotaions are bounded (and non-zero) at the crack tip has led to the 
formulation of new fracture criteria (Morozov, 1984). 

Due to the complexity of the equations of motion provided by the couple stress elastic theory with 
rotational inertia, only few crack propagation problems have been considered in the literature, most of them 
have been solved numerically and very few closed form solution have been worked out. In particular, Parhi 
& Das (1970) used the couple stress theory to investigate the stress distribution in a half-plane due to a 
moving line load on the surface. They assumed that the load is moving slower than both the longitudinal 
and transverse wave speeds of the elastic media (subsonic case) and found that the couple stress effect 
depends on the speed of the moving load, the couple-stress parameter, and the Poisson's ratio. Han et al. 
(1990) investigated dynamic propagation of a finite crack under mode-I loading in a micropolar elastic solid. 
Since the Wiener-Hopf technique can not be easily applied to finite crack problem, these authors solved 
numerically a pair of two-dimensional singular integral equations obtained by using an integral transform 
method. They provided solutions for dynamic stress intensity and couple stress intensity factors by using 
the obtained values of the strengths of the square root singularities in macrorotation and the gradient of 
microrotation at the crack tips. Itou (1981) considered the steady-state propagation in plane strain condition 
for the case of vanishing rotational inertia and only one characteristic length, analysing the influence of the 
crack tip speed on the asymptotics of the stress field. However, as already pointed out by Zhang et al. 
(1998) and Radi (2008) for the static case, a simple asymptotic characterization of the crack tip fields is not 
sufficient to analyse the fracture behaviour of materials with microstructure. This is due to the fact that the 
asymptotic solution is valid in a region close to the crack tip which is smaller than the characteristic lengths 
of the material, and thus it is of scarce physical relevance. A full field analysis is then required in order to 
grasp the qualitative and quantitative behaviour of the solution in a larger region and to be able to judge 
on the stress level supported by the material. 

The static full field solution obtained by Radi (2008) for a Mode HI crack by using Fourier transforms 
and Wiener-Hopf technique, shows that ahead of the crack tip within a zone smaller than the characteristic 
length in torsion, the total shear stress and reduced tractions occur with the opposite sign with respect to 
the classical LEFM solution, due to the relative rotation of the microstructural particles currently at the 
crack tip. However, this zone is found to have limited physical relevance and to become vanishing small 
for a characteristic length in torsion of zero. In this limit case, the solution recovers the classical A'ni field 
with square root stress singularity. Outside the zone where the total shear stress is negative, the full field 
solution exhibits a bounded maximum for the total shear stress ahead of the crack tip, whose magnitude 
was adopted as a measure of the critical stress level for crack advancing. The introduction of a stress based 
fracture criterion thus defines a critical stress intensity factor, which increases with the characteristic length 
in torsion. 

In this paper we analyse the steady-state propagation of a Mode HI crack in couple stress elastic materials 
with finite characteristic lengths in bending and torsion. The static analysis is extended to the case of steady- 
state propagation in order to study the effects of inertia and crack-tip speed on the stress and deformation 
fields, as well as the variation of the fracture toughness due to the presence of microstructure. The paper 
starts with the description of the fully dynamical version of the couple stress elastic model in Section 2, 
followed by an overview of some peculiar dynamical effects that this model is able to simulate. Differently 
from conventional elastic theories, this reasonably simple model can indeed predict the dispersive character 
of shear wave propagation through elastic media with microstructures, for a specific range of values of the 
microstructural parameters. Results provided in Section 2.1 concerning dispersive wave propagation can thus 
be used to validate the characteristic length scales by comparison with atomistic and molecular dynamics 
simulations (Maranganti & Sharma, 2007) or with experimental observations (Jakata & Every, 2008). Next, 



the problem of steady-state crack propagation under Mode III loading conditions applied to the crack surfaces 
is presented in Section 2.2. The analytical full-field solution is then addressed in Section 3, making use of 
Fourier transform and Wiener-Hopf technique and basically following the approach introduced by Atkinson 
& Leppington (1977) for a plane stationary crack problem. The crack tip is assumed to propagate at subsonic 
speed. Closed- form solutions are provided in Section 3.1 for the special case of vanishing rotational inertia, 
and in Section 3.2 for an arbitrary but small value of the rotational inertia parameter. 

The stability of antiplane crack propagation is discussed in Section 4 by using the fracture criterion 
based on the maximum shear stress hypothesis introduced by Radi (2008), which can find application in 
crack propagation also. It follows, indeed, that propagation is unstable if the maximum shear stress occurring 
ahead of the crack tip during steady propagation is increasing as the crack tip speed increases. 

2 Problem formulation 

Reference is made to a Cartesian coordinate system (0,a;i,X2, xs) centred at the crack-tip at time t = 0. 
Under antiplane shear deformation, the indeterminate theory of couple stress elasticity (Koiter, 1964) adopted 
here provides the following kinematical compatibility conditions between the out-of-plane displacement M3, 
rotation vector (p, strain tensor e and deformation curvature tensor x 

1 du3 1 du3 1 du3 1 9u3 

2 0x1 2 0x2 2 0x2 2 0x1 

Therefore, rotations are derived from displacements and the tensor field x turns out to be irrotational. 
According to the couple stress theory, the non-symmetric Cauchy stress tensor t can be decomposed into a 
symmetric part cr and a skew-symmetric part r, namely t = cr + t. In addition, the couple stress tensor 
fi is introduced as the work-conjugated quantity of x'^ ■ For the antiplane problem within the couple stress 
theory e, cr, x, x ^iid fi are purely deviatoric tensors. The reduced tractions vector p and couple stress 
tractions vector q are defined as 

p = t^n + -V/i„„ X n, q = fi'^n - /i,i„n, (3) 

respectively, where n denotes the outward unit normal and yu„„ ~ n ■ fjbn. The conditions of dynamic 
equilibrium of forces and moments, neglecting body forces and body couples, write 

9cri3 9(723 9ti3 (9r23 .. a//ii 9^21 , „ J.. dfii2 , 9/^22 „ J.. ,.. 

axi 0x2 oxi 0x2 oxi 0x2 oxi 0x2 

where p is the mass density and J is the rotational inertia. 

Within the context of small deformations theory, the total strain e and the deformation curvature x "^re 
connected to stress and couple stress through the following isotropic constitutive relations 

^ = 2Ge, ^Ji^2Ge{x'^lJx), (5) 

where G is the elastic shear modulus, I and ?/ the couple stress parameters introduced by Koiter (1964), with 
— 1 < r; < 1. Both material parameters I and rj depend on the microstructure and can be connected to the 
material characteristic lengths in bending and in torsion, namely 

4 = ^7^2, it = Wl + V- (6) 

Typical values of (,b and it for some classes of materials with microstructure can be found in Lakes 
(1986; 1995). The limit value of 77 = —1 corresponds to a vanishing characteristic length in torsion, which 



is typical of polycrystalline metals. Moreover, from the definitions ^ it follows that £t ~ it for 77 = —0.5 
and it = i = V^lb for rj = 0. The constitutive equations of the indeterminate couple stress theory do not 
define the skew-symmetric part r of the total stress tensor t, which instead is determined by the equilibrium 
equations (|1])2.3. Constitutive equations (O together with compatibility relations ([1]) and ^ give stresses 
and couple stresses in terms of the displacement u^: 
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Mil = -M22 = Gr(l + r;) 



Tl3 

where A denotes the Laplace operator. By means of ([7|) and ^, the equation of motion Qi becomes 

Gl? J 

GAU3 - -Y^^U3 + -Aug = pU3. (10) 

2.1 Preliminary analysis on shear wave propagation 

In this section we analyse the propagation of shear waves in a couple stress elastic material. We assume a 
solution of the governing equation ([TU|) in form of a planar shear wave as follows 



M3 



(a;i,a;2,i)=Ge'('="-"-'^*\ (11) 



where C is the amplitude, k is the wave number and ui the radian frequency. A substitution of pip into pop 
yields the following dispersive equation 



where Cg = yG/p is the shear wave speed for classical elastic material and 9 = \JAG/ J. By solving (|12p. 
we get the two equivalent dispersive relations 




+^r' "^f k^+i ' (13) 

where /i is a characteristic length given hy h = Cs/9. We obtain that the shear wave is dispersive and 
propagates with phase speed 
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2 k'^Phl + V 



(14) 



where h^ = h/i is a dimensionless parameter. Note that the relation p^ was found by Ottosen et al. (2000) 
for the special case of vanishing rotational inertia, /ig = 0. 

The normalized shear wave speed c/cs as a function of the normalized frequency uj(./ Cg is shown in Fig. [1] 
for different values of the parameter h^, which is a measure of the rotational inertia. In the case of vanishing 
rotational inertia, ha = 0, the wave propagation is dispersive with monotonically increasing phase speed, 



as found by Ottosen et al. (2000) (see also Askes & Aifaiitis, 2011). In the presence of rotational inertia, 
ho y^ 0, the dispersive character of the propagation is reduced until, for the limiting case ho = l/v2, the 
wave propagation is non-dispersive as in classical elastic materials. Beyond this limit, the wave propagation 
is again dispersive but with decreasing phase speed. It is worth noticing that the same result was found by 
Askes & Aifantis (2011) in the framework of stable gradient elasticity. Similar behaviour was also reported 
by Georgiadis & Velgaki (2003) for Rayleigh wave propagation in CS material. In particular, they analysed 
the behaviour of dispersive curves in the range 1/V24 < ho < -^/lO/S of the rotational inertia (or, in their 
notations, 0.25 < jj-h'^/fj < 20, where ft = G, fi = Gt"^ /2 and h^ = 3J/Ap). However, the authors of the 
aforementioned papers did not discuss the limit value for the parameter ho separating the two essentially 
different regimes. Nevertheless, analysing the numerical results reported by Georgiadis and Velgaki (Figs. 5 
and 6 in their paper), one could expect that the limit value coincides with our findings. Since this value is 
found for two different wave problems (antiplane shear waves in the present case and plane strain Rayleigh 
waves in their case), we believe that it may be a fundamental constant for CS materials. As it will be seen 
in Sec. 01 this limit value is also related to the stability of crack propagation. 




Figure 1 : Dispersive character of shear waves in couple-stress materials. 

It is also noted from Fig. [l]that, for low frequencies, cj — > 0, the shear wave propagates with the same 
speed as in classical elastic materials, c = Cs, because long wavelengths do not feel the microstructure. Vice 
versa, for high frequencies, w — > oo, short wavelengths interact with the microstructure and the shear wave 
propagates with constant speed depending on the rotational inertia, c = Cs/yjzho- Consequently, in the 
range of rotational inertia < /lo < l/v2, the crack propagation is subsonic provided that the crack tip 
speed V is less than Cg, whereas in the range ho > l/v^, the crack propagation is subsonic subject to the 
condition V < Cs/V^ho- 



2.2 Steady-state crack propagation 

We assume that the crack propagates with a constant velocity V straight along the xi-axis and is subjected 
to reduced force traction p3 applied on the crack faces, moving with the same velocity V, whereas reduced 
couple traction qi is assumed to be zero. 



P3ixi,0^,t)=TT{xi-Vt), qi{xi,0^,t)^0, for xi-Vt<0. 



(15) 



We also assume that the function r decays at infinity sufficiently fast and it is bounded at the crack tip. 
These requirements are the same requirements for tractions as in the classical theory of elasticity. 

It is convenient to introduce a moving framework A" = I'l — Vt, y = X2. By assuming that the out of 
plane displacement field has the form 

U3(xi,a;2,i) =ui(A, y), (16) 



then the equation of motion PH)) writes: 
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where rri = V/cg is the normahzed crack velocity. We consider the subsonic regime (see Fig. [5]), so that 
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Figure 2: Subsonic and supersonic regimes in the m — ho plane. 

According to ([3]) , the non- vanishing components of the reduced force traction and reduced couple traction 
vectors along the crack hne y = 0, where n = (0, ±1, 0), can be written as 



, , 1 ^^^22 
^3 = ^23 + -^^, gi=/i21, 



(19) 



respectively. By using ^2, (IH])i,2; ®2, and ([T9| . the loading conditions (fT5|) on the upper crack surface 
require the following conditions for the function w: 



dy 2 dy 
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Ahead of the crack tip, the skew-symmetry of the Mode III crack problem requires 



?7^;p^ = 0, for X<0, 2/ = 0+. (20) 



w = 0. 



dy^ 



V 



9X2 



0, for X >Q, y^O^ 



(21) 



Note that the ratio 77 enters the boundary conditions (pn)) -(pT |) . but it does not appear into the governing 
PDF dUl). 



3 Full field solution 

In the present section the Wiener-Hopf analytic continuation technique (Noble, 1958; Freund, 1990) is used 
to obtain the full field solution for a semi-infinite crack propagating in an infinite medium with constant 
velocity and subject to moving loading applied on the crack surfaces. Only the upper half-plane (y > 0) is 



considered due to the skew-symmetry of the problem. Use of the Fourier transform and inverse transform is 
made. For the function w(X, y) they are 

w{s, y) = p w{X, y)e*^^dX, w{X, v) = ^ J ^i^, y)e-''''ds, (22) 

respectively, where s is a real variable and the line of integration C will be defined later. 

Introduction of (P^ ? into the governing equation ([TT]), yields the following ODE for W{s,X2) 

- 1 [£^ (1 - 2n.^hl) .^ + 1] + ^ [f (1 - 2n.^hl) s^ + 2{l- m^) s^] W ^ 0. (23) 

Eq. (|23p admits as bounded solution in the upper half-plane 

W{s, y) = C(.s)e-"('^^)«/^ + D{s)e-P^'''^^y'\ for y > 0, (24) 

where 

a{z)^^l + {l-hlm^)z^ + x{z), P{z) = ^l + [l - hlm^)z^ - x{z). (25) 

X{z) = sjl + 2(1 - hl)m?z'^ + h^m^zJ. (26) 

The function xi^) has branch points at ±ibi and ±162, where 

In the case ho < l/v^, the branch points are located along the imaginary axis and the branch cuts are 
chosen to run from ibi to 162 and from —ibi to — i&2- In the opposite case, ho > l/v2, the branch points are 
symmetrically located in the four quadrants (separeted from the real axis) and the branch cuts are chosen 
to run from 161.2 to ioo and from ~ibi_2 to —ioo. In both cases the branch of the function x(^) is chosen 
such that the function is positive along the real axis. 

The function /3(z) has a branch point at z = 0, whereas the points ±ic, where 



2(1 -m2) , , 

are branch points of either a{z) or l3{z) depending on the values of the parameters m and ho- Consequently, 
in view of (fT8|) . c > and the additional branch cuts for the functions a{z) and /3(z) can be chosen to be 
along the imaginary axis. Again, the branches are chosen such that the functions are positive along the real 
axis. 

The functions C{s) and D{s) can be determined by imposing the boundary conditions ((20)) and pip . 
Applying the Fourier transform to the function w(X, 0) and keeping into account ()2ip i we obtain 

w{s,0)=W-{s), (29) 

which is analytic in the lower half complex s-plane, Im s < 0. The Fourier transform of the reduced force 
traction p^, and reduced couple traction qi at y ~ leads to 



dw e 



(2 + ,y-2m^/^g).^^(.,O) + 0(.,O) 



^ + ^, (30) 



o2 

^(s,0)+77s2w(s,0)=0, (31) 

where pg (s) is the Fourier transform of the reduced traction ahead of the crack tip and is analytic in the 
upper half complex s-plane, Ims > 0. 

A substitution of (p4|) into ((29|) . ([30| and ([3T|) allows us to find the functions C{s) and 15(5), in the form 



and also the following functional equation of the Wiener-Hopf type 



^(^) = ry— 52 ' (32) 



^(*) = -.2 _ m ' (33) 



pt{^) + Yi^{'^^)^-{s)^r-{s), (34) 

where 

f{z) = -L^laPia" + /32 + 27?z2) + ^2^2 _ ^2 ^i\ (35) 

a + p I. J 

In order to apply the Wiener-Hopf technique of analytic continuation, Eq. (j34p needs to be factorized into 
the product of two functions analytic in the upper and lower half planes, respectively. This is described, in 
the following sections, separately for the two cases of vanishing (J = 0) and non-vanishing (J 7^ 0) rotational 
inertia. 

3.1 The case of vanishing rotational inertia J = 

In the case of vanishig rotational inertia, we have J = and consequently h^ = 0. The Wiener-Hopf equation 
([M)l reduces to 

pt{s) + YiI{s(^)^-{s)=^r~{s), (36) 

where 

/(z) = ^-{a/3(a + /?)2- (a/3 -77^2)2}^ (37) 

a + p ^ 

and 

a{z) ^ Vl + z2 + x(^), /3(z) = ^Jl + z■^-x{z), x{z) = v^l + 2m2z2. (33) 

The function x{^) has branch points at ±i6, where h = l/{^/2ra), and the branch cuts are chosen to run 
from ±i& to ±icx) along the imaginary axis. 

The function /3(z) has a branch point at z = 0, whereas the points ±zc, where 



c = ^2(1-"^'), (39) 

are branch points of either a{z) or /3(z) depending on the values of the parameter m. Consequently, the 
additional branch cuts for the functions a{z) and P{z) can be chosen to be along the imaginary axis. 

It is noted that /(z) does not have any pole in the entire complex plane and has the following asymptotics 
at zero and infinity along the real line 

/(z) = i(l+7?)(3-77)|z|3 + 0(|z|2), z^ioo, zeR, (40) 



/(z) = 2v/l-m2|z| + 0(|z|2), z^O, zeK. (41) 

Moreover, the only zero of /(z) along the real axis is at z = 0. Note that 77 = —1 is a special case as 
the behaviour of the function /(z) at infinite differs and consequently the solution is different. We will not 
consider this limit case. 



3.1.1 Wiener Hopf factorization 

In order to apply the Wiener-Hopf technique of analytic continuation, eq. (j36p needs to be factorized into 
the product of two functions analytic in the upper and lower half-planes, respectively. 
Instead of the function f{z), we will factorize the function 

Hz) = -J^ (42) 

where 

*(z) = i(l + 77)(3 - T])z^ + 2V1 - m2. (43) 

It is noted that, along the real axis, the function k{z) has no zeros and has the following asyniptotics at 
zero and infinity 

fc(z) = l + 0(|z|), z^O, zeR. (44) 

fc(z) = l + 0(|z|-^), z->±oo, zeR, (45) 

Moreover, the function k{z) has two simple poles at the points ±id, where 



d = 2^^^^. (46) 

Vi^ + V)i3-V) 

The function k{z) is factorized as follows 

fc(^)-^, zeR, (47) 

where fc"*"(z) and k~{z) are analytic in the upper and lower half-planes, respectively, and are defined as 
follows 

e^^'^ for lmz>0, 

fc+(z) = { ^Biz) (48) 

— , for lmz = 0, 

e«(^) for lniz<0, 

k-{z)={ ^ ^^ (49) 

v/fc(I)e-"(^' for lmz = 0, 



R{z) = 


2m J_^ t - z 


TTl Jo t^ - Z2 


z 

■Kl 


\j:t->^ j:t:i>] 



(50) 

where use is made of the condition fc(— z) ~ k{z). The integrals in ((50|) must be evaluated as the Cauchy 
principal value if Im z = 0. 

The functions fc+(z) and k^ (z) have the following asyniptotics at zero and infinity along the real axis 

fc±(z) = l-— -/ ^ ^ ' dt + 0{z^), z^O, zeR, (51) 



2Tri 1^ t2 



1 /"^ 

fc±(z)==H / logk{t)dt + 0{z~^), z^zboo, zeM. (52) 

27riz J_^ 



In the next section, we will also need to factorize the function v z^. This is accomplished by writing 

y^^z^J^/'^y', (53) 

where the branch cut for z^ ranges from to =pioo. The functions z^ and z_ take along the real line 
the limit values 



,1/2 



1/2 
^+ = 



z>0, 



1/2 



i[-zfl\ z<0, 



zl/2, Z > 0, 

-z(-z)l/2, z<0. 



(54) 



3.1.2 Solution of the Wiener Hopf equation 

By making use of P^ . the Wiener-Hopf equation (P5|) becomes 

G 



p+(s) + y^2£2[(l + ,^)(3 _ ^)s2^2 _^ ^^jY~n^-\-k{st)W~ {s) = T-(s). 



The term \/ s^P is factorized as follows 



^f7^^{sllf_l'^{sifl'^, 



(55) 



(56) 



where the functions (s^)^ and (s^)_ are analytic in the upper and lower half-plane and are defined in Sec. 



By using (|47| and ((56|). eq. ([55| can be written as 



^^^^^^ + S(^^)^/^[(l + .)(3 - r,)sH- + AVT^]k-isi)W-i.s) . Ilif)^ 



We assume the following form of the loading applied on the crack faces 



is£)X- 



r{^X) = ^e'"\ X<0, 



that is. after Fourier transform, 



so that eq. ([57]) becomes 



T-{S) 



Tn 



1 + isL' 



T^k+{si) 



[si)^ *<- \s")+ (l + *sL) 

The right-hand side of (|60p can be easily split in the sum of plus and minus functions, namely 



nk+{s(.) 



Tn 



{si)]^\l + isL) l + tsL 



k+{s£) 



1 + isL ' 



where 



„ k+{U/L) 



[zi/L) 



1/2- 



(57) 

(58) 
(59) 

(60) 



(61) 



(62) 



The first term of the right-hand side is no longer singular at s = i/ L, and thus it is a plus function, whereas 
the last term is a minus function. 



10 



By using (|FT1) and (|T7)) . eq. ([50)1 can be factorized as follows 



fc+(s£)p3+(s) Tc 







(s^)i/2 1 + isL 



^+ 



k+{se) ^ 



C / „xl/2r/-, , \/o X 2^2 , . A 7i,-/ /,\ / N , ^0'^ 



Asi)]/^[{1 + ?y)(3 - T])s^f + 4\/l - m^]k- {s£)w'' (s) + - . 



(63) 



The left and right hand sides of (|63p are analytic functions in the upper and lower half-planes, respectively, 
and thus define an entire function on the s-plane. The Fourier transform of the reduced force traction ahead 
of the crack tip and the crack opening gives p^ ^ s^ " and w^ ^ s^^'"^ as s — > oo. 

Therefore both sides of ()63p are bounded as s — > oo and must equal a constant F in the entire s-plane, 
according to the Liouville's theorem. 

As a result, we obtain 

-+(^ ^0 ^„(f^l^£(l_MsL) 

— ( ^ 2roS£ l-F{\ + isL) 

w [s) = -—^ 775 , (65) 



where the function ^ is defined in (j43|) . 

The constant F is determined by the condition that the displacement w{X) is zero at the crack tip X — Q, 
that is 



/ 

J —I 



w-{s)ds = 0, (66) 

J —00 

which leads to 

ds 



^ = — ^ dS • (^^) 



-00 {s(.f_!^'^{se.)k-{se.) 

3.1.3 Analytical representation of displacements, stresses and couple-stresses 

The reduced force traction ahead of the crack tip P3{X) and the crack opening w{X) can be obtained from 
(|64| and ([65]) . respectively, by inverse Fourier transform, according to ([22|l 9. Since the integrand does not 
have branch cuts along the real line, the path of integration C coincides with the real s-axis. Further, we 
introduce the change of variable ^ = s£, thus obtaining 



The Fourier transform of stress (symmetric and skew-symmetric) and couple stress fields can be derived from 
dMl), dSH), ([331) and dZl)-®, namely 

0-23(5,0) = -2To^ —■ YR, ' ('0) 

" + ^ {sefr{\^isV)^{si)k-{si) 
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T23(s,0) = -ToS^-f (1 + 7y)/£4 + ^2^2 [^2 ^ 2r; + (1 + 7^)^^^^/sH^ + c-']]: 
a + n L J 



a + /3 
l-F(l + isL) 



{se.fl^{\ + isL)^{sl)k-{siy 



7l22(s,0) = -2roS£(l + 77)M) 



Q,^ - ,ys2£2 



I- F{l+isL) 



The inverse Fourier transform can be performed as explained above, thus obtaining for X > Q 



aMO) = -^ 



<om - ve i~F{i+ i^L/e) 



ixiji 



dt 



-— ■ /"OO 



r23(^,0) 



Ai22(X,0) 



'27r^ 



{TqE 






«(C)+/5(0 

1 - F(l + i^i/^) 



eJ\i+iwmm-io 



(i + v) 



e 






'ixiji 



di. 



(71) 



(72) 



(73) 



(74) 



(75) 



3.1.4 Results 



The variation of the normalized total shear stress ahead of the crack tip at y = 0, namely t23^/2o, versus 
normalized distance Xji is shown in Fig. |3l for three different values of r; = {—0.9, 0, 0.9} and three different 
values of the normalized crack tip speed m = {0.01,0.5,0.99}. The steady-state solution approaches the 
static solution for small values of 771 (see the solid lines in Fig. [3] corresponding to 771 — 0.01), confirming 
earlier results obtained by Radi (2008) for a stationary crack. In particular, the present results confirm that 
the strain-gradient effect is confined in a region of size hi close to the crack tip, where the solution gradually 
changes from the classical positive field, acting at large distance, to the negative near-tip asymptotic field. 
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Figure 3: Variation of total shear stress t23 along the X-axis, ahead of the crack tip. 

As discussed in Radi (2008), the zone where the shear stress has the negative sign is very small and can 
be considered of no physical importance. Outside this zone, the total shear stress attains a maximum, t^23^^ 
which is bounded and positive for — 1 < 7/ < 1 and < 777 < 1. As the crack tip speed increases, the total 
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shear stress tends to decrease ahead of the crack tip, and the rate of decrease is more pronounced for positive 
and large values of the ratio 77, namely for large values of the characteristic length in torsion according to 
(p])?. Correspondingly, the zone ahead of the crack tip where the total shear has the negative sign enlarges 
and may reach the size of 2€, where ^ is proportional to the characteristic length in bending according to 

®i. 

The corresponding normalized variations of the symmetric stress, skew-symmetric stress and couple stress 
are shown in Figs. 21in|and|n|for the same values of r] and m. In Fig. |5|it can be observed that the symmetric 
shear stress 1723 is positive and finite ahead of the crack tip and tends to display the singular behavior of 
classical elasticity as 77 tends to —1. A different trend is observed as the crack tip speed becomes faster for 
negative or positive values of 77. Specifically, the symmetric stress increases with the crack speed for negative 
values of 77, whereas it decreases for large and positive values of 77. 
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Figure 4: Variation of symmetric shear stress a2i along the X-axis, ahead of the crack tip. 

The skew-symmetric shear stress T23 ahead of the crack tip is negative and strongly singular (Fig. O. 
Its magnitude increases with the crack tip speed for every value of 77, so that the total shear stress always 
decreases as the crack tip speed becomes faster, as already observed. 
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Figure 5: Variation of skew-symmetric shear stress r23 along the X-axis, ahead of the crack tip. 

The plots in Fig. |6| show that the couple stress 7.^22 decreases with X and tends to vanish for X ^ i, 
thus recovering the classical solution of linear elasticity. The decrease is very fast for negative values of the 
ratio rj. For example, for 77 — —0.99 then 7^22 completely vanishes at X = 5£, whereas for positive values of 
T], namely for larger value of the characteristic length in torsion, the effects of couple stress and rotational 
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gradients are felt at larger distance from the crack tip. In particular, the non-monotonic behaviour of fj.22 
already observed in the static analysis (Radi, 2008) is recovered for very slow crack tip speed, but it tends 
to disappear for faster speed of propagation. 
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Figure 6: Variation of the couple stress /i22 along the X-axis, ahead of the crack tip. 

The normalized variations of the sliding displacement on the crack face, wG/Tq at y = 0, versus normal- 
ized distance X/£ to the crack tip are shown in Fig. |7| for the selected range of values of 77 and m. As well 
known, the crack tip profile is blunted for the classical elastic solution, but it turns out to be sharp for CS 
elastic materials. The magnitude of the sliding displacement between the crack faces slightly decreases as rj 
increases from —1 to 1, thus indicating that the crack becomes stiffer due to the effects of the microstructure, 
here brought into play by strain rotational gradients. This trend is confirmed also for crack propagation and 
proves that the microstructure may shield the crack tip from fracture, as already observed for stationary 
crack (Zhang et al., 1998; Georgiadis, 2003; Radi, 2008). The results plotted in Fig. [Tjalso show that the 
magnitude of the crack opening profile increases as the crack speed becomes faster, for every set of material 
parameters. Since the crack tip sliding displacement is expected to be proportional to the square-root of 
the generalized J-integral (see Piccolroaz et al., 2012), then an increasing amount of energy is required in 
order to increase the crack tip speed, especially as the crack tip speed approaches the shear wave speed, thus 
predicting stable crack propagation within the full subsonic regime. 



3.0 

2.5 


--^ 


1] = -0.9 




;, =0 




1] = 0.9 

m = 0.01 


2.0 
1.5 








~"v^ 




--- m =0.8 
""-.^ m =0.99 


1.0 
0.5 
0.0 


--- 


'^'''^ 




''''^^^^^ 




^^:::::-:X^ 


5 


-4 -3 -2 -1 C 


) - 


5 -4 -3 -2 -1 C 


) - 


5 -4 -3 -2 -1 C 






X / ( 




X 1 ( 




X 1 ( 



Figure 7: Variation of the crack opening displacement to along the crack face. 

The results of the full field analysis for very small values of X plotted in Figs. |5H7| agree with the 
asymptotic solution obtained by Piccolroaz et al. (2012). However, the radius of validity of the asymptotic 
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solutions becomes vanishing small as 77 approaches —1, namely as £t — > 0. In this limit case, the classical 
elastic solution is indeed recovered, since both the symmetric and the total shear stresses become unbounded 
at the crack tip and display the classical square-root singularity. 

The variations of the non dimensional ratios ^23^/21), X'^^'^^/l and X'^ /£ with 77 are plotted in Fig. [5] 
There, it can be observed that the maximum shear stress decreases as rj increases, as well as the crack tip 
speed becomes faster. Correspondingly, both the distance to the crack tip of the location of the maximum 
shear stress and the size of the zone ahead of the crack tip with negative shear stress become larger. 
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Figure 8: Variations with the ratio r\ of the maximum shear stress t^^^ and its location X™^'' and size X" of the zone with 
negative shear stress ahead of the crack tip. 



3.2 The case of non- vanishing rotational inertia J 7^ 

In the case J 7^ 0, the asymptotics of the function /(z) at zero and infinity along the real axis are as follows 

j{z)^T{^Mm)V? + 0{\zf), z->±oo, zeM, (76) 

(77) 



/(z) = 2v/l-m2|z| + 0(|zn, z^O, zeK, 



where 



T{ri,hom) = 



l-jf - 2h^m^ + 2^1 -2hlm^{l + tj - hlm^) 



(78) 



1 + y/l - 2hlm'^ 

Note that for 77 = and any horn < l/-\/2 this constant is always positive. However, if 77 7^ 0, there exists 
a set of parameters in the plane (77, horn) where T(77, ham) < for all hgrn < 1/V2, see Fig. [9l As for the 
case of vanishing rotational inertia J = and 77 = —1 described above, we are not considering the case when 
T(77,/7om)<0. 

It is convenient to factorize the function k{z) in the form (j42p where ^(z) is different and defined as 
follows: 

*(z) = T(r7, hom)z^ + 2^1 - m^- (79) 

One can show that, under the condition T(ri,hom) > 0, the procedure for factorization is analogous to 
that reported in the previous section, and also the expressions of displacements, stresses and couple-stresses 
(except the skew-symmetric stress) can be represented in the same form as in Sec. 13.1.31 
The Fourier transform of the skew-symmetric stress takes the form 

T23(s,0) = -ToE^—-\a^l3^ + {a^ + /3^ + aP)'ns^f - (1 - 2hlm^)s^£^{i-is^£^ - a/3)]x 
a + p I i 

l-F(l + isL) 

^ Tn ^ ' 

(si) ' {1 + isL)'ii{sl)k- (si) 



(80) 
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which, after inversion, gives 



o^{i)+m) 



{«'(0/3'(e) + (a'(0 + P\o + a{i)m))rie 



(1 - 2him')e{ve - a{om))]^w- 



1 - F(l + i^L/f.) 



^-iXill 



cr{\^iiLii)nm-{i) 



dt 



(81) 



3.2.1 Results 



The effect of the rotational inertia on the stress, couple stress and displacement fields ahead of the crack 
tip can be observed in Figs. [T0Hl4l There, the normalized variation of these fields are plotted for the same 
crack tip speed, namely for m — 0.8, for three different values of ry = {—0.9, 0, 0.9} and three different values 
of the non-dimensional parameter ho = h/i providing for the contribution from rotational inertia. It can be 
observed that the effects of rotational inertia on these fields mainly consist in the rising of the stress level 
ahead of the crack tip as the rotational inertia parameter ho is increased. 

However, for large value of rj the symmetric shear stress display peculiar behavior just ahead of the crack 
tip, where it decreases and becomes negative for ho > l/v^- Correspondingly, the couple stress /i22 also 
change its sign ahead of the crack tip. As it will be discussed in the next section, such behavior occurs as a 
consequence of the loss of stability of the crack propagation. 

Moreover, Fig. [14] shows that the crack opening profile slightly increases with ho thus denoting that an 
additional amount of energy must be supplied from the external loadings to allow the crack propagating 
at constant speed if the rotational inertia contribution is considered. Correspondingly, the maximum shear 
stress increases its magnitude and tends to concentrate nearer to the crack tip, as illustrated in Fig. 1151 
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Figure 10: Variation of total shear stress i23 along the X-axis, ahead of the crack tip. 




Figure 11: Variation of symmetric shear stress (T23 along the X-axis, ahead of the crack tip. 
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Figure 12: Variation of skew-symmetric shear stress T23 along the X-axis, ahead of the crack tip. 
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Figure 13: Variation of the couple stress /^22 along the X-axis, ahead of the crack tip. 



1.0 
0.8 

t? 0.6 

O 

^ 0.4 

0.2 

0.0 



■=e.c. '' " 


-0.9 & m 


= 0.8 : 




N;<\ 


\ 


ho 


= 0.01 


\. 


— ho 


= 0.2 


V 


-— ho 


= 0.28 


\: 



-3 -2 
X I ( 





'7 


= & m 


= 0.8 




'"""--T^* 


""^s. ^^ 






- 










- 






^ 




- 




ho 


= 0.01 


\\ 




— 


ho 


= 0.707 


\\ 




— 


ho 


= 0.8 




\ 



-3 -2 
X I ( 



'1 


= 0.9 &m 


= 0.8 


~--C~-""-:; 






- 


\^ 


\\ 


ho 


= 0.01 


Vo, 


ho 


= 0.707 




-■- ho 


= 0.8 


\' 



-3 -2 
X I ( 



Figure 14: Variation of the crack opening displacement w along the crack face. 
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4 Stability analysis of the crack propagation 

In this section we discuss the stabihty of the crack propagation on the basis of the full field solution presented 
above. It must be pointed out that a sound criterion of fracture propagation in materials with microstructures 
cannot be based on the classical parameters of the crack-tip stress singularities. The singular behaviour occurs 
indeed at a distances to the crack tip smaller than the characteristic lengths, which can be commensurate 
with the microstructural size and even smaller. This would contradict the basic concept of continuum 
modelling which cannot predict phenomena occurring at the scales smaller than the microstructural size. 
According to that, Dyskin et al. (2010) recently proposed an intermediate asymptotic analysis for a finite 
length crack corresponding to the distances from the crack tip larger than the characteristic internal lengths 
but yet smaller than the crack length. 

In the present work, we adopt the fracture criterion of the maximum total shear stress, introduced in 
Radi (2008) and Geogiadis (2003), which usually occurs at distance commensurate with the microstructural 
size and, thus, still in the field of validity of continuum mechanics. The occurrence of a maximum positive 
value of the total shear stress ahead of the crack tip, ifjf , allows us to formulate a simple criterion by 
assuming a critical shear stress level tc at which the crack starts propagating. The fracture criterion is then 
formulated as follows 

t^r = rc, (82) 

where tc is a material parameter. 

The normalized variations of the maximum attained by the total shear stress t™^^, its location ahead of 
the crack tip and the size of the zone with negative shear stress versus the normalized crack tip speed m 
are plotted in Figs. \W[ [T71 and [T51 for rj = {—0.9,0,0.9}, respectively. The results here reported show that 
for values of the rotational inertia parameter ho larger than l/-\/2 the maximum shear stress tf^^ grows 
very rapidly and becomes unbounded at a crack-tip speed lower than the shear wave speed Cs, namely for 
771 < 1. This behaviour necessarily implies a limit speed of propagation lower than Cs and also unstable 
crack propagation if the stress based fracture criteria (|82|) is assumed. Conversely, if the rotational inertia 
contribution is small, namely for Hq ^ 1/ v2, and rj is sufficiently greater than —1, then the maximum shear 
stress decreases as the crack speed becomes faster and tends to zero as the crack speed approaches Cs (solid 
lines in Figs. ll7H18p . Therefore, within this range of values for rj and ho, the crack propagation turns out to 
be stable, according to the fracture criteria (|82t . Non-monotonic behaviour is observed in the plot of tgjf 
versus m for rj approaching the limit value —1 (solid lines in Fig. 1161 plotted for rj ~ —0.9). In this case, 
the maximum shear stress first decrease as the crack tip speed increases, thus suggesting the occurrence of 
stable crack propagation at speed sufficiently lower than the shear wave speed. However, as m approaches 1 
then the curve rapidly increases. 

Therefore, the presence of microstructures provides a stabilizing effect on the propagation of crack. 
Indeed, for large values of the characteristic length in torsion, namely for large positive values of rj, cracks 
can propagate steadily at high speed, also for large contribution of the rotational inertia. Differently, for 
small characteristic length in torsion, namely for negative values of eta, crack propagation becomes unstable 
as soon as the contribution of the rotational inertia becomes more pronounced. In particular, for vanishing 
small contribution from rotational inertia, crack propagation seems to be stable also for small size of the 
microstructures. 

These findings in part agree with the results provided by Kennedy & Kim (1993). These authors per- 
formed finite element analysis for the case of a crack propagating at a constant velocity in a micropolar 
elastic material and observed that the micropolar effects consist in reducing the energy release rate for low 
velocity cracks and became less pronounced with increasing crack propagation velocity. 
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Figure 16: Variations with the crack speed m of the maximum shear stress f^jf" ^'^^ ^^^ location X™'"' and size X" of the 
zone with negative shear stress ahead of the crack tip. 
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Figure 17: Variations with the crack speed m of the maximum shear stress i^''" and its location X™''" and size X of the 
zone with negative shear stress ahead of the crack tip. 
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Figure 18: Variations with the crack speed m of the maximum shear stress i^''" and its location X™'"' and size X of the 
zone with negative shear stress ahead of the crack tip. 



20 



5 Conclusion 

The full field solution for a semi-infinite crack steadily propagating in a couple stress elastic material and 
subject to antiplane loading on the crack surfaces has been obtained by using Fourier transform and the 
Wiener-Hopf technique. The material microstructure is characterized by finite characteristic lengths in 
bending and torsion. The effect of rotational inertia is also considered in the analysis. 

The main conclusion is that couple stress elasticity, providing for material characteristic lengths, may 
provide a clearer picture of the failure process near a rapidly moving crack tip than classical elasticity. The 
analysis based on classical theory dictates that the only way a Mode III crack can respond to additional 
influx of energy is by accelerating until it eventually outruns the energy at the shear wave speed. However, 
phenomena such as the actual limiting crack speed under Mode III loading condition, which has been observed 
to be much lower than the shear wave speed as predicted by classical elastic theory, cannot at all be addressed 
by classical, linear elastic theories. Experimental investigations (Livne et al., 2010) also bring to light the 
influence of the fracture process on the limiting value. Since the process zone near the crack tip is strongly 
influenced by the microstructural parameters, such as the material characteristic lengths and the rotational 
inertia, the maximum crack speed is also expected to be influenced by them. The fact that a fracture process 
zone has structure and dynamics associated with its evolution provides for the dependence on microstructural 
parameters for the fracture energy and presents a possibility for explaining the lower observed limiting speeds 
and the definition of the range of stable crack propagation. Recently, advanced models of the fracture process 
have been proposed in literature, ranging from idealized atomistic and lattice models to phenomenological 
models to mechanism-based nucleation and growth models. While none of these models have advanced to the 
stage of providing quantitative explanation of the limiting speed or for the dependence of the fracture energy 
on the crack speed, there are qualitative features of the present micropolar model that appear promising. 
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